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Amendments to the Specification : 

Please replace paragraph [001 1] as follows: 

[001 1] In order to simulate a well test, whatever the medium, this equation has to be 
solved in space and in time. D e f i n i ng -Dlscretization o f the reservoir (mesh pattern) is 
therefore performed and solution of the problem fieds -consists in finding t ha pressures 
of the meshes with time, itself defieed -discretized in a certain number of time intervals. 

Please replace paragraph [0019] as follows: 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0019] Other features and advantages of the method according to the invention will be 
clearfrom reading the description hereafter of a non limitative realisation example, with 
reference to the accompanying drawings wherein : 

- Fig. 1 shows an example of a 2D fracture mesh, E M Fracture M esh. 

- Fig. 2 shows an example of a 2D matrix block, M S Matrix Block . 

- Fig. 3 shows an imposed flow rate variation curve F(t) in a well -test, 

- Fig. 4 shows an example of a fracture network for which the method according to the 
invention leads to a substantial reduction in the number of meshes to be processed, 
and 

- Figs. 5 and 6 illustrate flow charts of the method of the present invention. 

" — Fiqs 7 an ^ 8 illustrate the input data as describe?- in the Assignee's U.S. Patent 
6.023.656. 
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Please replace paragraph [0041] as follows: 
[0041] Apart from the pressures, the terms of this equation are known. The pore 
volumes of the fracture meshes and of the matrix blocks (ft) are known by means of 
the mesh pattern, the fracture-fracture and matrix-fracture transmissivrtres (Tg) are 
calculated as described above and flow rates (Q ; ) are zero everywhere except at the 
well-reservoir connections where they are imposed. The mosh pattefFMS-kftewB 
using the method d e scribed in th e abovo citod US p a tent 6 t 02 3 t 6 lli 6, Fracture 
le ngths - Q Bd Op e n i ngs (widths) nrn giwan hy t h n r h ?ra c to ristics of t he fra c t uro 
network and arc input data> - T he pore volume in a given fracture mesh is determined 
from a volumetric calculation using th«* length -dimensions of fractures within the 
mesh-ao44 ho correspond i ng op e ning of th e fracturo . 

Please replace paragraph [0057] as follows: 
[0057] The geometry of the fracture network and the attributes of the fractures 
(conductivity, opening) are given in the form of a file as in the method described in 
the aforementioned US patent 6,023,656. 

[0057a1 As shown on FIGS. 7 and 8 fractures are assumed to be { substantially 
vertical (i.e. perpendicular to the iaver limitslHowever. a same data structure can 
applied to fractures slightly deviati ng from the vertical direction. The 3D image is 
discretized vertically complying with th e actual geological layering if such inform^ nn 
i§, available. If not, any arbitrary dis cretization is applied to the Imag * F af ^ 
horizontal layer L Is characterized h y its vertical coordinate 7L in th* > referent 
system of coordinates (OX. OY. Q 7). 
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r0057b1 For each layer L. a series of rectangles R has fobe defin e F a ch rectangle 
consists in a fracture plane elem ent comprised between the limits , of a given lave r. 
Hence, each natural fracture consists in a se t of superimposed rec ta ngles R and is 
assigned an origin (fracture origin). Each rectangle is defined bv: 

the three coordinates (xC). vO. zO) of the rectangle origin Q, For a given 
natural fracture, all the origin points of the constitutive rectang les are situated on th e | 
same vertical (or highest dip) line drawn from the fracture orig in; 

the co-o rdinates of the horizontal unit vector i (xH. V H) and nf th ^ vertical unit 
vector j fxV, vV) defining the orient ation of the rectangle in the refe renc e system of 
co-ordinates, with x Vertical and v Ve rtical being zero in case of vertical fractures hut 
considered as input d ata to be able to deal with non-vertical fractu res; 

the two algebraic horizontal lengths 1 - and 1+ separating the origin of the 
rectangle and the tw o lateral (vertical) limits of this rectangle: 

theheiqht h of the rectangle, that is the length of the renton ^ ^n Vn 
direction I which is the laver thic kness if discretization along direction i fits the 
geology: 

the hydraulic conductivity c d erived from the application of H arry* l«w to 
fracture flow (for a pressure gradient 

/ 
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the How rate in the fracture with a height h Is 

ch_ AP 

4_beinq the fluid viscosity ). The conductivity c is given bv the relation c-ka 
where k=a ; /12 fusing PoiseniNe's Idealized representation of fractures) is the intrinsic 
permeability of the fracture and a its eouivaient hydraulic aperture a. The hydraulic 
conductivity c is a ref erence value given for a direction of the maximum horizontal 
stress para l lel to the fracture direction; 

the two uppe r and lower neighboring rectangles UR. LR: 

the fracture set to which the rectangle considered belongs to ; 

theo rjentation angle a ? of the direction of maximum horizontal stress taken 
from (OX) axis in the reference system of coordinates: 

for ench fracture set, a correlation table correlating 1 ) the ar>gfe between the 
direction of maximum horizontal stress and fracture direction fazimuth) with 2) the 
hydraulic qonductivitv c or equivalent hydraulic aperture a previous l y defined. 
"Horizontal" and "vertical" stand in the context for directions respectively parallel and 
perpendicular to the limits of layer s which hare are assumed horizontal. Layer limits 
discrettze fracture planes in the 'Vertic a l" direction. It must be pointed out that the 
aforesaid i nput data 1) are suitable for all the existing software tools used for 
characterizing and generating fra c ture and 2) could be used to fl>nretlze a network 
of slightly non-vertica l fractures, i.e. not perpendicular to layer limit? , 
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